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' ABSTRACT: * 

This paper reports investigation on laminar free convection in a vertical double passage channel for 
electrically conducting fluid in the presence of first order chemical reaction. The channel is divided into 
two passages by inserting a thin plane conducting baffle. After placing the baffle one of the passage is 
concentrated. An analytical solution has been developed for the coupled nonlinear ordinary differential 
equations using regular perturbation method. The results show that the thermal and mass Grashof 
number and Brinkman number enhances the flow where as Hartmann number and first order chemical 
reaction parameter suppresses the flow at all the baffle positions in both the streams. 
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I. INTRODUCTION 

The study of heat and mass transfer for an electrically conducting fluid under the influence of 
transverse magnetic field has attracted researchers in the recent past due to its relevance in many engineering 
problems such as MHD generators, heat exchangers, thermal processing, gas-cooled nuclear reactors and others. 
In the absence of magnetic field, references [1-3] will give some ideas about fluid flow and thermal 
characteristics inside a vertical channel with symmetric or asymmetric thermal boundary conditions. For an 
infinite vertical plate, Raptis and Kafoussias [4] studied the flow and heat transfer characteristics in the presence 
of magnetic field. Later Raptis [5] extended the vertical plate problem to a vertical channel problem in the 
presence of magnetic field. Malashetty et al. [6, 7, 8], Umavathi et al. [9], and Prathap Kumar et al., [10] studied 
mixed convection in a vertical channel for one and two fluid models. Umavathi et al. [11] analyzed magneto 
hydrodynamic free convection flow in a vertical rectangular duct for laminar, fully developed regime taking into 
consideration the effect of Ohmic heating and viscous dissipation. Umavathi and Sridhar [12] studied the 
Hartmann two-fluid Poiseuille-Couette flow in an inclined channel. Shah and London [13] have summarized the 
laminar forced convection heat transfer results for various channel cross sections, which were reported in the 
literature until 1970. 

In particular process involving the mass transfer effects has been considered to be important precisely in chemical 
engineering equipments. The other applications include solidification of binary alloys and crystal growth dispersion of 
dissolved materials or particulate water in flows, drying and dehydration operations in chemical and food processing plants, 
evaporation at the surface of water body. The order of the chemical reaction depends on several factors. One of the simplest 
chemical reactions is the first order reaction in which rate of reaction is directly proportional to the species concentration. 
Das et al. [14] have studied the effect of mass transfer on the flow started impulsively past an infinite vertical plate in the 
presence of wall heat flux and chemical reaction. Muthucumaraswamy and Ganeshan [15, 16] have studied the impulsive 
motion of a vertical plate with heat flux/mass flux/suction and diffusion of chemically reactive species. Seddeek [17] has 
studied the finite element method for the effect of chemical reaction, variable viscosity, thermophoresis, and heat 
generation/absorption on a boundary layer hydro magnetic flow with heat and mass transfer over a heat surface. Kandasamy 
et al. [18, 19] have examined the effects of chemical reaction, heat and mass transfer with or without MHD flow with heat 
source/suction. The rate of heat transfer in a vertical channel could be enhanced by using special inserts. These inserts can be 
specially designed to increase the included angle between the velocity vector and the temperature gradient vector, rather than 
to promote turbulence. This increases the rate of heat transfer without a considerable drop in the pressure by Guo et al. [20]. 
A plane baffle may be used as an insert to enhance the rate of heat transfer in the channel. To avoid a considerable increase 
in the transverse thermal resistance into the channel, a thin and perfectly conductive baffle is used. The effect of such baffle 
on the laminar fully developed combined convection in a vertical channel with different uniform wall temperatures has been 
studied analytically by Salah El-Din [21]. Their study showed that for mixed convection the heat transfer between the walls 
and fluid can be significantly enhanced according to the baffle position and higher values of Nusselt number can be obtained 
when the baffle become as near the wall as possible. 
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Keeping in view the applications of chemical reaction and the increase of rate of heat transfer by 
introducing a baffle, motivated to investigate the effect of first order chemical reaction for electrically 
conducting fluid in a vertical channel. After inserting the baffle the fluid in stream-I is concentrated. 

II. MATHEMATICAL FORMULATION 

Consider a steady, two-dimensional laminar fully developed free convection flow in an open ended 
vertical channel filled with purely viscous conducting fluid. The X-axis is taken vertically upward, and parallel 
to the direction of buoyancy, and the F-axis is normal to it as seen in Fig. 1 . The channel walls are maintained at 
a constant temperature and the fluid properties are assumed to be constant. The channel is divided into two 
passages by means of thin, perfectly conducting plane baffle and each stream will have its own pressure gradient 
and hence the velocity will be individual in each stream. 



Stream-I 



Stream-II 



Y = -h 



Y = h 



Figure 1. Physical configuration. 
The governing equations for velocity, temperature and concentrations are 



Stream-I 

p g J3 T ( 7i - Tw 2 ) + p g J3 C ( C - C w 2 ) - ^ + M - a e Ft] U , = 0 



+ — I ~ I + — B: U 

dY 2 k\dY J k 

D^-KC = 0 
dY 2 

Stream-II 

P g/3 T (T 2 -T W2 )-- + M — 



, a e B„U 2 =0 



d% jLif dU 2 V ct. 



0 



dY 2 ' k{dY ) k 

Subject to the boundary and interface conditions on velocity, temperature and concentration as 



(1) 

(2) 
(3) 

(4) 
(5) 



U x =0, T i= T W] 
U 2 =0, T, 



C = C W1 at Y = -h 
at Y = h 



AT AT 

U,=0, U 2 =0, T,=T 2 , —L = ^-, C = C W1 at Y = h" 
2 1 2 dY dY 



(6) 
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Introducing the following non-dimensional variables, 



U, „ T,-T W2 

u = = , a, = — 

II T -T 

u 1 ± W I 1 IV, 



Gr- 



gjB T ATh 3 



Gc 



gj3 c ACh 3 



u~ 



u~ 



C -C 
C -C 



„ U,h 
Re = — , Br 
v 



kAT ' h 



h dp . „ „ „ v , a B n h 
£_ AT = T -T AC = C -C Y-— M =— - - 

juU l dX h ju 



(7) 



One obtains the non-dimensional momentum, energy and concentration equations corresponding to 
stream-I and stream-II as 



Stream-I 

d 2 u, 



+ GR T 9 X +GR c <j)-p-Mu l =0 



- + Br 



dy 2 
d 2 e, 

~dy 

— T -a <p = 0 

dy 2 



'duA 2 , , 
— - I +M u; 
dy 



(8) 
(9) 
(10) 



Stream-II 

d 2 u 



dy 2 

d 2 0 2 



+ GR T 0 2 -p-M 2 u 2 =0 



- + Br 



du 2 

dy 



-M 2 u 2 



= 0 



Subject to the boundary conditions, 

t/j = 0 , 6> = 1 , (j) = 1 , at y = -1 
u 2 = 0 , 9 2 = 0 , at y = 1 

<f<9. d0, 

«j = 0 , u 2 = 0 , 0 x -0 2 , — L = — - , q> = n , at j = j*, 

dy dy 

where G/? r = — ,GR C = — ,a= , n=— —. 

Re Re D C 1 -C w2 



(11) 
(12) 



(13) 



III. SOLUTIONS 

Solution of equation (10) can be obtained directly using boundary condition (13) and is given by 



B l Cosh(ay} + B 2 Sinh(ay} 



(14) 



Equations (8), (9), (11) and (12) are coupled non-linear differential equations. Approximate solutions 
can be found by using the regular perturbation method. The perturbation parameter Br is usually small and 
hence regular perturbation method can be strongly justified. Adopting this technique, solutions for velocity, 
temperature and concentration are assumed in the form 



{y) = u io {y) + Bru il (y) + Br 2 u i2 (y) + ... 
0, {y) = e i0 {y) + Br0 il (y) + Br 2 9 i2 (y) + ... 



(15) 
(16) 



Substituting equations (15) and (16) in equations (8), (9), (11) and (12) and equating the coefficients of 
like power of Br to zero and one, we obtain the zero and first order equations as 
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Stream-I 

Zeroth order equations 



d 2 u 



dy 
d 2 0 



+ GR T 9 m +GR c (j)-p-Mu l0 =0 



ML 

2 ull r "10 1 ""c v f m ™-m 

k = o 



dy 

First order equations 

d 2 u x 
~dy 
d 2 0 l 
~df 



2 +G/? r 6> u -M k, 



du l0 

dy 



+ M 2 u u 2 



(17) 
(18) 

(19) 
(20) 



Stream-II 

Zeroth order equations 

d ^ + GR T e 2Q -p-M 2 k 20 =0 

dy 

First order equations 



d 2 u 2 

d 2 0 2l 
dy 2 



+ GR T 0 2l -p-M 2 u 2l =0 



+ M 2 m 20 2 



= 0 



(21) 
(22) 

(23) 
(24) 



The corresponding boundary conditions reduces to 
Zeroth-order 



-1 



d9 lQ d0 2Q 
0 2n , —— = ——, 0 = n, at y = y* 
dy dy 



u m = 0 , 0 1O = 1 , 0 = 1 , at y : 
«2o =0, 0 20 = 0, at y = 1 

u l0 = 0 , m 20 = 0 , 6 10 — ^ 20 

First order 

Mjj =0, 6> n =0 at y = -1 
m 21 = 0 , ^j, =0 at y = \ 

d0 df) 

u n =0, u 21 =O,0 n =0 2l , -7^ = -^ aty = y* 

dy dy 



(25) 



(26) 



The solutions of zeroth and first order equations (17) to (24) using the boundary conditions as in 
equations (25) and (26) are 
Zeroth order 
Stream-I 

0 W =z l y+z 2 



u l0 = A, Cosh{My) + A, Sinh(My) + r x + r 2 y + r 3 Cosh(ay) + r 4 Sinh(ay) 

Stream-II 

0 2O =z 3 y + z 4 

u 20 = A 3 Cosh(My) + A 4 Sinh(My) + r 5 +r 6 y 



(27) 
(28) 

(29) 
(30) 
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First order 
Stream-I 

6 U = E 2 + E t y + q x y 2 + q 2 y 3 + q 3 y 4 + q 4 Cosh{ay) + q 5 Sinh(ay) + q 6 Cosh(2ay) + q 1 Sinh(2ay) 

+q B Cosh(2My) + q 9 Sinh(2My) + q l0 Cosh(My) + q u Sinh(My) + q n yCosh(My) + q n ySinh{My) 
+q u yCosh(ay) + q l5 ySinh(ay) + q lft Cosh (« + M ) y + q xl Cosh [a -M ) y + q lg Sinh (or + M ) y (31) 
+q l9 Sinh(a - M}y 

u n = E 5 Cosh(My) + E 6 Sinh(My) + H l + H 2 y + H 3 y 2 +H 4 y 3 +H 5 y 4 +H 6 Cosh{ay) + H 1 Sinh{ay) 

+H i Cosh(2ay) + H 9 Sinh(2ay) + H l0 Cosh(2My) + H n Sinh(2My) + H l2 yCosh(My) + H l3 ySinh(My) 
+H u y cosh(ay) + H l5 ySinh(ay) + H l6 Cosh [a + M ) y + H n Cosh (a-M)y + H ls Sinh [a + M ) y (32) 
+H l9 Sinh(a -M)y + H 20 y 2 Cosh(My) + H 2l y 2 Sinh(My) 

Stream-II 

0 2l =E A + E 3 y + Fj 2 + F 2 y 3 + F 3 y 4 + F 4 Cosh(2My) + F 5 Sinh(2My) + F 6 Cosh(My) + F^inhiMy) 

+ F s yCosh(My) + F 9 ySinh(My) 
u 2i = E 7 Cosh(My) + £ 8 sinh(My) + H 22 + H 23 y + H 24 y 2 + H 25 y* + H 26 y 4 + H 21 Cosh(2My) + H 2S Sinh(2My) 

+ H 29 yCosh(My) + H 30 ySinh(My) + H 3l y 2 Cosh(My) + H 32 y 2 Sinh(My) 

The constants appeared in the above equations are presented in the Appendix. 

IV. RESULTS AND DISCUSSIONS 

The problem of free convective heat and mass transfer in a vertical double passage channel filled with 
electrically conducting fluid is investigated. The analytical solutions are found using regular perturbation 
method considering Brinkman number as the perturbation parameter. The effects of governing parameters such 
as Hartmann number, thermal Grashof number, mass Grashof number, Brinkman number and first order 
chemical reaction parameter on the velocity, temperature and concentration are shown graphically. 

The effect of Hartmann number on the flow is shown in Figs. 2a,b,c and 3a,b,c. As the Hartmann 
number increases, the velocity and temperature decreases in both the streams at all the baffle positions. As the 
Hartmann number increases the fluid decreases which is the classical Hartmann result. When the baffle is near 
the hot wall the maximum velocity is in stream-II, when the baffle is near the cold wall the maximum velocity is 
seen in stream-I and when the baffle is in the center of the channel the maximum velocity is seen in stream-I. 

The effect of thermal Grashof number on the velocity and temperature field is shown in Figs. 4a,b,c 
and 5a,b,c at all three different baffle positions. As the thermal Grashof number increases, the velocity and 
temperature increases in both the streams at all the baffle positions. This is an expected result because increase 
in thermal Grashof number results in increase of bouncy force and hence increases the flow in both the streams 
at all the baffle positions. 

The effect of mass Grashof number on the velocity and temperature field is shown in Figs. 6a,b,c and 
7a,b,c respectively. As the mass Grashof number increases the flow is enhanced in both the streams at all the 
baffle positions. The enhancement on velocity is significant in stream-I when compared to stream-II .This is due 
to the fact that the fluid is concentrated only in stream-I. Though the fluid is not concentrated in stream-II still 
its effect is observed in stream-II when the baffle position is in the center of the channel, this is due to the reason 
that we have considered the baffle to be conducting. That is say that, there is heat transfer from stream-I to 
stream-II but there is no mass transfer from stream-I to stream-II. Due to transfer of heat from stream-I to 
stream-II results in increase of both thermal buoyancy force and concentration buoyancy force and hence 
velocity increase in stream-II slightly as mass Grashof number increases as seen in Fig. 6. As mass Grashof 
number increases the temperature increases significantly in stream-I when the baffle position is at the center of 
the cold wall. There is no much variation on the temperature field when the baffle position is near the hot wall. 

The effect of Brinkman number is shown in Figs. 8a,b,c and 9a,b,c on the velocity and temperature 
fields respectively. As the Brinkman number increase both the velocity and temperature increases in both the 
streams at all the baffle positions. This is due to the fact that increase in Brinkman number increases the viscous 
dissipation and hence the flow is enhanced. 



(33) 
(34) 
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The effect of first order chemical reaction parameter on the velocity, temperature and concentration 
fields are displayed in Figs. 10a,b,c, lla,b,c and 12a,b,c respectively. As a increases the velocity, temperature 
and concentration decreases in stream-I, and remains constant in stream-II. The similar result was also obtained 
by Srinivas and Muturajan [22] for mixed convective flow in a vertical channel. This is due to the fact that the 
fluid in stream-I is concentrated. The maximum value of velocity and temperature is seen in stream-II for the 
baffle position at y* = -0.8 and in stream-I for the baffle position at y* = 0 and 0.8. 

V. CONCLUSIONS 

The effect of first order chemical reaction in a vertical double passage channel filled with electrically 
conducting fluid by inserting a thin baffle is investigated. The following conclusions were drawn: 

1. Increasing the values of Hartmann number reduces the flow field where as increase in the thermal 
Grashof number and mass Grashof number enhance the flow in both the streams at different baffle 
positions. 

2. Increase in the perturbation parameter (Brinkman number) enhances the velocity and temperature in 
both the streams. 

3. Increase in the chemical reaction parameter suppresses the velocity, temperature and concentration in 
stream-I and remains invariant in stream-II. 

REFERENCES 

[1] S. Habchi and Archarya. Laminar mixed convection in a symmetrically or asymmetrically heated vertical channel. Numer. Heat 
Transfer, 9:605-618, 1986. 

[2] W. Aung and G. Worku. Developing flow and reversal in a vertical channel with asymmetric wall temperatures. J. Heat Transfer, 
108:299-304, 1986. 

[3] M. Iqbal, B.D. Aggarwala and A.G. Fowlere. Laminar combined free and forced convection in a vertical non-circular ducts under 

uniform heat flux. Int. J. Heat Mass Transfer, 12:1123-1139, 1969. 
[4] N. Raptis, Kafoussias. Heat transfer in a flow through a porous medium bounded by an infinite vertical plane under the action of 

magnetic field, Energy Res., 6:241-245, 1982. 
[5] Raptis. Flow through a porous medium in the presence of magnetic field. Energy Res., 10:97-100, 1986. 

[6] M.S. Malashetty, J.C. Umavathi and J. Prathap Kumar. Two fluid magnetoconvection flow in an inclined channel. Int. J. 

Transport Phenomena, 3:73-84, 2000. 
[7] M.S. Malashetty, J.C. Umavathi and J. Prathap Kumar. Convective magnetohydrodynamic two fluid flow and heat transfer in an 

inclined channel. Heat and Mass Transfer, 37:259-264, 2001. 
[8] M.S. Malashetty, J.C. Umavathi and J. Prathap Kumar. Magnetoconvection of two immiscible fluids in a vertical enclosure, Heat 

and Mass Transfer, 42:977-993, 2006 
[9] J.C. Umavathi, B. Mallikarjun Patil and L Pop. On laminar mixed convection flow in a vertical porous stratum with symmetric 
wall heating conditions. Int. J. Transport Phenomena, 8:1-14, 2006. 
[10] J. Prathap Kumar, J.C. Umavathi and Basavaraj M Biradar. Mixed convection of composite porous medium in a vertical channel 

with asymmetric wall heating conditions. J. Porous Medium, 13:271-285, 2010. 
[11] J.C. Umavathi, LC, Liu, J. Prathap Kumar and I. Pop. Fully developed magneto convection flow in a vertical rectangular duct. 

Heat and Mass Transfer, 47 : 1 - 1 1 , 20 1 1 . 
[12] J.C. Umavathi and K.S.R. Sridhr. Hartmann two fluid Poisulle-Couette flow in an inclined channel, Int. J. Appl. Mech. Engg., 
14:539-557, 2009. 

[13] R.K. Shah and A. London. Laminar flow forced convection heat transfer and flow friction in straight and curved ducts-A 

summary of analytical solutions. Stanford University, Stanford, California, USA, Tech. Rep. vol. 75, 1971. 
[14] U.N. Das, R.K. Deka and V.M. Soundalgekar. Effects of mass transfer on the flow past an impulsively started infinite vertical 

plate with constant heat flux and chemical reaction. Forsch Ingenieurwes, 60:284-287, 1994. 
[15] R. Muthucumaraswamy and P. Ganesan. On impulsive motion of a vertical plate with heat flux and diffusion of chemically 

reactive species. Forsch Ingenieurwes, 66:17-23, 2000. 
[16] R. Muthucumaraswamy and P. Ganesan. First order chemical reaction on flow past impulsively stratified vertical plate with 

uniform heat and mass flux. Acta Mech., 147:45-57, 2001. 
[17] M.A. Seddek. Finite element method for the effects of chemical reaction, variable viscosity, thermophoresis and heat 

generation/absorption on a boundary layer hydro magnetic flow with heat and mass transfer over a heat surface. Acta Mech., 

177:1-18,2005. 

[18] R. Kandasamy, K. Perisamy and K.K. Sivagnana Prabhu. Effects of chemical reaction, heat and mass transfer along a wedge with 
heat source and concentration in the presence of suction or injection. Int. J. Heat Mass Transfer, 48:1388-1394, 2005. 

[19] R. Kandasamy, K. Perisamy and K.K. Sivagnana Prabhu. Chemical reaction, heat and mass transfer on MHD flow over a vertical 
stretching surface with heat source and thermal stratification effects. Int. J. Heat Mass Transfer, 48:4557-4561, 2005. 

[20] Z.Y. Guo, D.Y. Li and B.X. Wang. A novel concept for convective heat transfer enhancement, hit. J. Heat Mass 
Transfer,.41:2221-2225, 1998. 

[21] M.M. Salah El-Din. Fully developed laminar convection in a vertical double-passage channel. Appl Energy, 47:69-75, 1994. 
[22] S. Srinivas, R. Muthuraj. Effect of chemical reaction and space porosity on MHD mixed convective flow in a vertical asymmetric 
channel with peristalsis. Mathematical and Computer Modeling, 54:1213-1227, 2011. 



www.ijceronline.com Open Access Journal Page 6 



Effect of First Order Chemical Reaction on Free Convection in . 



NOMENCLATURE 





magnetic field 


Br 


Brinkman number 


c, 


concentration in Stream-I 


c 0 


reference concentration 


c P 


specific heat at constant pressure 


c 

p 


dimensionless specific heat at constant pressure 


Q 


Ulll UalVJll V_<V_JG 1 llv_ 1G IIIl) 


En 


applied electric field 


E 


electric field load parameter 


h 


channel width 


h* 


width of passage 


8 


acceleration due to gravity 


Gr 


Grashoff number 


Gr 


Grashoff number 


G 


modified Grashoff Number 


GR T ,GR C dimensionless parameters 


k 


thermal conductivity of fluid 


M 


Hartmann number 


P 


nondimensional pressure gradient 


Re 


Reynolds number 


T„T 2 


dimensional temperature distributions 


T ,T 


temperatures of the boundaries 




reference velocity 




nondimensional velocities in Stream-I, Stream-II 


u v u 2 


dimensional velocity distributions 


y* 


baffle position 


GREEK SYMBOLS 


a 


chemical reaction parameters 




coefficients of thermal expansion 


Pc 


coefficients of concentration expansion 


a e 


electrical conductivity 


AT, AC 


difference in temperatures & concentration 




non-dimensional temperature 


7 


kinematics viscosity 


$ 


non-dimensional concentrations 


P 


density 


M 


viscosity 



SUBSCRIPTS 

i refer quantities for the fluids in Stream-I and Stream-II, respectively. 
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Figure5 .Temperature profiles for different values of thermal Grashof 
number GR at (a) \ ' =-(). 8 (to) y*=0 (o) y*=0.8 
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Figure6.Velocity profiles for different values of concentration 
Grashof number GR, at (a) y*=-O.S (to) v*=0 (c) y*=0.8 
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Figure 1 1 .Temperature profiles for different values of chemical 
reaction parameter a (a) y*=-0.8 (b) y*=0 (c) y*=0.8 



www.ijceronline.com Open Access Journal Page 12 



Effect of First Order Chemical Reaction on Free Convection in . 




www.ijceronline.com Open Access Journal Page 13 



Effect of First Order Chemical Reaction on Free Convection in . 



Appendix 

Sinh(ay) + nSinh(a) nCosh(ay ) — Cosh(ay ) 1 

£?] = ; ; — , B 2 = ; — , Zj = 

Cosh(a)Sinh(ay ) + Cosh(ay )Sinh(ay ) Cosh(a)Sinh(ay ) + Cosh(ay )Sinh(ay ) 2 

_j_ -1 _ GR T C 2 -P _ GR T C I GR^ GR C B 2 _ GR T C 4 -P 

z 2 --, z 3 --, Zt- 2 , 1 - M 2 . ^ M * ' ^""^T^' ^ " M 2 ' 

GR, u 3 * 
r 6 = M 2 ' ' ^ =-ri+r 2 -r 3 cosh(a) + r 4 sinh(a), r 3 =_ r5 J 4 =-r 5 _ rg j , 

T.Sinh(My) + T 2 Sinh(M) 

T- = —(r. + r 9 y + n coshfay ) + r, sinh(ay )) , A = , 

Cosh{M)Sinh(My) + Cosh(My)Sinh(M) 

T 2 Cosh{M)-T 1 Cosh{My) T 3 Sinh(My*)-T 4 Sinh(M) 
^2 = 7^ , 777777; , 777 * N 777, 7777777777 77777 ' ^3 ~~ 



Cosh{M)Sinh{My ) + Cosh{My )Sinh{M) ' 3 Cosh{M)Sinh(My ) - Cosh(My )Sinh(M) ' 

= _3 3 ^ V A / =-2MV 2 , p 3 =-M 2 r 2 , 

4 S&ift(flf) 2 F2 12 3 2 

/a \ 2 ('"3 2 +'-4 2 )(« 2 +^ 2 ) , 2 2 , 

p 4 =-(2M r 1 r 3 +2r 2 r 4 aj, p 5 =2M r { r 4 +2r 2 r 3 a, p 6 =-- -> P 7 =-(M r 3 r 4 + r 3 r 4 a J, 



p s = -(A 2 M 2 +A 2 2 M 2 ), p 9 =-2A l A 2 M 2 , p l0 =-(2A l M 2 r l +2A 2 Mr 2 ), p u = -(2A 2 M 2 r 1 +2A l Mr 2 ), 
p 12 = -2\r 2 M 2 , p l3 = -2A 2 M 2 r 2 , p 14 = -2M 2 r 2 r 3 , p 15 = -2M 2 r 2 r 4 , 

p 16 = -(\r 3 M 2 + A 2 Mr 4 a + A 2 M 2 r 4 + A l Mar 3 ), p„ = -(A i M 2 r 3 + A 2 Mar 4 - A 2 M 2 r 4 -A l Mar 3 ), 



Pi 



p is = -(A 2 r 3 M 2 + A i Mr 4 a + \M 2 r 4 + A 2 Mar 3 ) , p 19 = -(A i M 2 r 4 + A 2 Mar 3 -A 2 M 2 r 3 -A i Mar 4 ), 9i = ~ 
p 2 _p 3 _p 4 a-2p l5 _p 5 a-2p l4 _ p 6 _ p 7 



l2= — > ft =777' ft = 5 , q 5 = - , q 6 =—j, q 7 =—5-. ft =7777 < ft 



6 12 14 « 3 " « 3 " 4« 2 ' 4cT 4M 2 4M 

_ p 10 M-2p l3 _ p u M-2p l2 = p }1 _ = p ll _ _p u _p j± _ p lb 

fto ..3 ' ftl ,,3 ' ft2 ,,2 ' ft3 ..2 > ftt _ 2 ' fts _ 2 ' ftf> / M ,\2 

M M M M a a (a + M) 



fts = , P ' 8 7j , = , Pl9 , 2 . ^,=-(M 2 r 5 2 +M 2 r () 2 ), /? 2 =-2MV6. » 3 =-M\ 2 , 
(a-M) (a+M)' (a-M) v ; 

R 4 =-A 2 M 2 -A 4 2 M 2 , R 5 =-2A 3 A 4 M 2 , R 6 = -(2A 3 A/V 5 +2A 4 M 2 r 6 ) , 7? 7 = -(2A 4 M 2 r 5 + 2A 3 M 2 r 6 ) , 



7?L F _7?l F = 7?1 F _ ^4 p _ _^5_ p _ R 6 M-2R 9 _ R 7 M-2R & _R^ = R^_ 
2' 2 ~6 3 12' 4 4M 2 ' 5 4M 2 ' 6 M 3 ' 7 M 3 ' 8 M 2 M 2 



f 



q i -q 2 + q 3 + q 4 Cosh(a) - q 5 Sinh(a) + q b Cosh(2a) - q 1 Sinh(2a) + q g Cosh(2M ) - q 9 Sinh(2M)^ 
+q w Cosh(M) — q n Sinh(M) — q l2 Cosh(M) + q l3 Sinh(M ) — q l4 Cosh(a) +q 15 Sinh(a) 
+q ib Cosh(a + M J + q l7 Cosh (a-M)- q lg Sinh (a+M)- q 19 Sinh (a-M) 



T 6 = -(F l +F 2 +F 3 + F 4 Cosh(2M) + F 5 Sinh(2M) + F 6 Cosh(M) + F 7 Sinh(M) + F s Cosh(M) + F 9 Sinh(M)) , 

T 7 = F^y" 2 + F 2 y" 3 + F 3 y" 4 + F 4 cosh(2My*) + F 5 Sinh{2My") + F b Cosh{My") + F 7 Sinh{My*) + F g y*Cosh(My*) 
+F 9 y" Sinh{my *) - (q l y* 2 + q 2 y*^ + q 3 y* 4 + q 4 Cosh{ay ) + q 5 Sinh{ay ) + q b Cosh(2ay ) + q 7 Sinh{2ay*) 
+q % Cosh{2My , ) + q 9 Sinh(2My*) + q w Cosh{My") + q n Sinh(My *) + q^'CoshiMy") + q u ySinh{My*) 
+q 14 yCosh(ay ) + q l5 ySinh(ay ) + q l6 Cosh(a + M)y + q 17 Cosh(a -M)y 
+q K Sinh(a + M)y* + q i9 Sinh(a - M)y") 

T 8 = 2Fj + 3F 2 y* 2 + 4F 3 y' 3 + 2MF 4 Sinh(2My) + 2MF 5 Cosh(2My') + F b MSinh(My) + F 7 MCosh{My) 
+F^(y*MSinh(My*) + Cosh(My )) + F 9 (y" MCosh(My") + Sinh(My')) - (2q l y* + 3q 2 y" 2 + 4q 3 y 3 
+q 4 aSinh(ay*) + q 5 aCosh(ay ) + 2aq b Sinh{2ay*) + 2q 7 aCosh{2ay*) + 2Mq s Sinh(2My*) 
+2Mq 9 Cosh(2My) + q 10 MSinh(My*) + q n MCosh(My ) + q l2 (y MSinh(My ") + Cosh(My")) 
+q 13 (y' MCosh(My*) + Sinh(My*)) + q u (y aSinh(ay*) + Cosh(ay")) + q l5 (y"aCosh(ay*) 
+Sinh(ay*))+ q ]b (a + M)Sinh(a + M)y* + q l7 (a - M)Sinh(a - M)y* + q u (a + M)Cosh(a + M)y* 
+q 19 (a- M)Cosh(a -M)y ) 
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_ (-T 5 +T 6 +T 7 +T s -T s y) _ (T 5 +T 6 + T 7 +T, -T s y) _ (-T 5 +T 6 + T 7 -T s -T g y') 
2 2 2 

(r s +r 6 -r 7 +r 8 +r 8 /) gr t (e 2 m 4 +2q { M 2 +2Aq 3 ) gr t (e i m 2 +6q 2 ) 

E A — 1 H, = , = , 

2 1 M M 

GR T (q 1 M 2 +l2q 3 ) GR T q-, „ Gfl 7 4 3 „ 2aq l5 GR T -q 4 GR T (a 2 -M 1 ) 

H 3 = ; , W , — ; — , H s — ; — , //, = z , 

M 4 M M (a 2 -M 2 f 

_ 2q 14 aGR T -g 5 GR T (a 2 -M 2 ) g 6 G/? r g 7 Gfl r „ _<ifi^T_ „ _%GR^ 

( a i-M 2 f 8 4« 2 -M 2 ' 9 4a 2 -M 2 ' 10 3M 2 ' 11 3M 2 

_ g n GR T — 2q n GR T M _ q l3 GR T -2q w GR r M q^GR, q^GR, 
n n — rm > "ii rr~i > n u = ; ;r , W,« - 



AM 2 ' 13 4M 2 ' 14 « 2 -M 2 ' 15 « 2 -M 2 ' 

g l6 GR T q^GRj q n GR T q, 9 GR T 

16 / s 2 t ' 17 , ,t -, ' 18 . ,7 , ' "19 , ,2 

(a + M)'-M 2 (a-M)'-M 2 (a + M)'-M 2 (a-M) -M" 

g l3 GR T q l2 GR T _ GR t (E 4 M 4 +2F 1 M 2 +2AF 3 ) _ GR T (E 3 M 2 +6F 2 ) 

H20 ~ AM ' 21 " AM ' ^ 22 " M 6 ' " 23 ~ M 4 

_ GR T (F l M 2 +12F 3 ) = GR I F ± = GR I F ± GR T F 4 GR T F 5 

M 4 ' 25 M 2 ' 26 M 2 ' " 3M 2 ' 28 " 3M 2 ' 

F S GR T -2F 1 GR T M F 9 GR r -2F 6 GR r M F 9 GR T F S GR T 

tl-, „ — ; , /I™ — r , tl„ — , tl. 



AM 2 ' 30 4M 2 ' 31 4M ' 32 AM ' 

T 9 =H l -H 2 +H 3 -H 4 + H s + H 6 Cosh(a) - H 7 Sinh{a) + H s Cosh(2a) - H 9 Sinh{2a) + Z/ 10 cosh(2M) 
-H u Sinh(2M) - H l2 Cosh{M) + H n Sinh(M) - H l4 Cosh{a) + H l5 Sinh(a) + H l6 Cosh(a + M) 
+H ll Cosh(a-M)- H n Sinh(a + M)- H l9 Sinh(a-M) + H 2Q Cosh(M) - H 2l Sinh(M) 

T 10 = H l +H 2 y' +H 3 y" 2 + H 4 y' 3 + H 5 y' 4 + H 6 Cosh(ay) + H 7 Sinh(ay* ) + H 9 Cosh(2ay) + H 9 Sinh(2ay) 
+H u) Cosh(2My) + H n Sinh(2My) + H \ 2 y 'Cosh(My ') + H l3 y*Sinh(My *) + H i4 y*Cosh(ay) 
+H ls y*Sinh(ay*) + H 16 Cosh(a + M)y* +H 17 Cosh(a-M)y* +H m Sinh(a + M)y* + H 19 Sinh(a-M)y* ' 
+H 20 y* 2 Sinh(My*) + H 2i y* 2 Sinh(My*) 

r u = H 22 +H 23 y* +H 24 y 2 + H 25 y" 3 + H 26 y' 4 + H 27 Cosh{2My) + H 2S Sinh(2My) + H 29 yCosh{My) 
+H 30 y*Sinh(My) + H 3l y 2 Cosh{My) + H 32 y 2 Sinh{My) 

T XJ = H 22 + H J3 + H 24 + H 25 + // 26 + H 27 Cosh(2M) + H 7S Sinh(2M) + H 29 Cosh(M) + H 30 Sinh(M) 
+H 3l Cosh(M) + H 32 Sinh(M) ' 

(T 9 Sinh{My ) + T IQ Sinh(M)) (T m Cosh{M)- T 9 Cosh{My )) 

E 5 = ~ „ ~ „ » E 6 ~ 



E, 



Cosh(M )Sinh(My ) + Cosh(My )Sinh(M) Cosh(M )Sinh(My ) + Cosh{My )Sinh(M ) 

(r i2 Sinh{My )-7j, Sinh(M ) ) fa, Cosh(M ) - T 12 Cosh(My* ) ) 

7 Cosh(M)Sinh{My )-Cosh{My )Sinh{M) ' 8 Cosh(M)Sinh(My)-Cosh(My)Sinh(M) 
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